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The existence of phantom energy in a universe which evolves to eventually show a big rip doomsday 
is a possibility which is not excluded by present observational constraints. In this letter it is argued 
that the field theory associated with a simple quintessence model is compatible with a field definition 
which is interpretable in terms of a rank-three axionic tensor field, whenever we consider a perfect- 
fluid equation of state that corresponds to the phantom energy regime. Explicit expressions for the 
axionic field and its potential, both in terms of an imaginary scalar field, are derived which show 
that these quantities both diverge at the big rip, and that the onset of phantom-energy dominance 
must take place just at present. 

PACS numbers: 14.80.Mz, 98.80. Cq 



I. INTRODUCTION 

Cosmology has become the current commonplace 
where the greatest problems of physics are being concen- 
trated or clearest manifest themselves. Along with the 
nature of dark matter, the origin of supermassive black 
holes and several unexplained huge amounts of released 
astrophysical energy are only some illustrative examples. 
However, perhaps the biggest problem of all of physics 
be the so-called dark energy problem [1] by which it is 
known that nearly the seventy percent of the total energy 
in the universe is in the form of an unobserved vacuum 
energy which is responsible for the present accelerating 
expansion of the universe [2,3] and whose nature remains 
being a mystery today [4] . There has been a rather huge 
influx of papers in recent years trying to shed some light 
onto the existence and the kind of the possible stuff that 
may make up dark energy [5]. 

Four main candidates to represent dark energy have 
been hitherto suggested: A positive cosmological con- 
stant [6], the quintessence fields (which may [7] or may 
not [8] be tracked) , some generalizations of the Chaply- 
gin gas [9] , and the so-called tachyon model of Padman- 
abhan et al. [10]. Even though such models could be 
all accommodated to present observational constraints, 
these models are posing new or traditional problems in 
such a way that none of them become completely sat- 
isfactory. In particular, cosmological data from current 
observations do not exclude [11] but may be suggesting 
[12] values of the parameter u> in the perfect-fluid equa- 
tion of state p = up of the most popular quintessence 
models which are smaller than -1. If this were the case, 
then dark energy would become what is now named phan- 
tom energy [13], an entity violating the dominant energy 
condition and thereby allowing the natural occurrence of 
wormholes and ringholes and even their corresponding 
time machines in the universe [14]. 

Caldwell, Kamionkowski and Weinberg have recently 
noticed [15] that in the framework of quintessence models 
phantom energy may lead to a doomsday for the universe 
- which would take place at a big rip singularity - once 



clusters, galaxies, stars, planets and, ultimately, nucleons 
and leptons in it are all ripped apart. Even though such 
a big rip does not take place in some phantom-energy 
models that use Chaplygin-like equations of state [16], 
it seems to be the simplest and most natural possibility 
stemming from phantom energy. On the other hand, the 
big rip cosmic scenario interestingly adds an extra quali- 
tative feature to the known set of cosmological models as 
it introduces a curvature singularity, other than that at 
the big bang, at a finite, nonzero value of the cosmic time. 
However, in spite of the feature that the blowing up of 
the scale factor appears to be unavoidable in models with 
equation of state p = ujp and u> < —1, a complete account 
of the nature and properties of the field theory associated 
with such models has not been done yet. This paper aims 
at investigating the characteristics of the massless scalar 
field which allow emergence of a big rip singularity in 
the simplest of such phantom energy models, assuming 
a perfect-fluid equation of state. It will be shown that 
the stuff making up phantom energy can be interpreted 
to be a vacuum sea of cosmic axions which can be de- 
scribed in terms of the kind of rank-three tensor field 
strengths predicted in supergravity and strings theories. 

This paper can be outlined as follows. In Sec. II we 
argue in favor of the idea that superlight axions are the 
source of cosmic phantom energy. A simple cosmological 
model accounting for a big rip singularity in the case 
that vacuum is filled with phantom energy is discussed 
in Sec. III. Sec. IV contains the solution of the phantom 
field theory within the simple cosmological model of Sec. 
III. We check that both the scalar field and its potential 
also have a singularity at the big rip. The results are 
summarized in Sec. V. 



II. AXIONS AS THE SOURCE OF PHANTOM 
ENERGY 



The definition of the massless scalar field cf> which is as- 
sumed to make up dark energy is usually taken to be the 
conventional simplest one; that is, in terms of pressure p 
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and energy density p, 

P=\tf + v{4>), v = l -^-v^\ (l) 

where V(4>) is the field potential. From the equation of 
state p = top and Eq. (1) it immediately follows that 



being the scale factor), is added; i.e. in the flat geometry 
case 
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Thus, if the weak energy condition p > is taken to be 
always satisfied [17], the requirement p + p < 0, ui < —1 
from phantom energy in this kind of models necessarily 
implies that the field <j> ought to be pure imaginary; that 
is to say, if the vacuum energy density for a phantom vac- 
uum as referred to any timclikc observer has to be posi- 
tive, then the massless scalar field making up the phan- 
tom stuff should be pure imaginary. I will argue in what 
follows that, in the classical framework, such a massless, 
pure imaginary scalar field actually represents an axion 
describable as a rank-three antisymmetric tensor field, 
considering after the associated cosmic theory for such a 
field. In fact, the Lorentzian action that couples such an 
axion field to gravity can generally be written as 



■L r 



(3) 



where L m is the Lagrangian for observable matter and 
A 2 = A llva A llva if we choose the axion to be given by a 
three-form A = dB field strength so that dA = 0, that is 
as a rank-three antisymmetric tensor field strength of a 
type arising in supergravity-theory motivated quantum- 
gravity solutions [18]. The equations of motion derived 
from action (3) are 
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in which A 2 — A^ a pA^ the asterisk denotes Hodge 

dual, and is the momentum-energy tensor for ordi- 
nary observable matter. One now can check that any 
explicit solution to these equations of motion subject 
to the usual Friedmann-Robertson- Walker (FRW) sym- 
metry for the metric and the corresponding spherically 
symmetric ansatz for the axion antisymmetric tensor 
A = f(r)e (with e the volume form which, when inte- 
grated over a surface of constant radius, yields the area 
of the unit three-sphere A = 2ir 2 /T(2)), if we set e.g. 
L m = and f(r) = Const., is the same as the solution 
obtained from the equations of motion derived from the 
FRW action integral describing coupling of a real mass- 
less scalar field <f> to gravity where an extra boundary 
term accounting for the axion properties, — a 3 </></>|o (» 



where ./V is the lapse function. The key point then is that 
the solution to the above equations of motion is even also 
preserved when we omit such an extra boundary term 
in the action for the scalar field cj>, provided this scalar 
field is simultaneously rotated to the pure imaginary axis, 
<j> — ► i<&. If the scalar field is equipped with a potential 

V{4>), then an extra term — J Q T dtNA s V{<j>) should be 
added to the action. In that case, the axionic action 
can be obtained by simply rotating <j> — > i$ in an action 
containing the above potential term, without any extra 
boundary term. 

Such an property, which was first noticed in Eu- 
clideanized solutions such as wormholes and other in- 
stantons describing nuclcation of baby universes [19], is 
actually independent of the metric signature and, there- 
fore, applies also to our Lorentzian cosmological context. 
Thus, one can interpret that the stuff that makes up 
phantom energy can be regarded to be a rank-three an- 
tisymmetric tensor axionic field. It is usually thought 
that axions can explain the absence of an electrical dipolc 
moment for the neutron and thereby solve the so-called 
strong CP problem [20]. The axions are chargeless and 
spinlcss particles with very tiny mass which interact with 
ordinary matter only very weakly. Such particles are be- 
lieved to have been abundantly produced in the big bang. 
It is worth noticing that whereas relic axions are an excel- 
lent candidate for the dark matter in the universe [21], 
their vacuum quantum background could make cosmic 
phantom energy. 

The coincidence and fine tuning problems could be 
thought to become exacerbated in the present scenario 
where one sets u> constant and < — 1. However, in dark 
energy models such as the generalized Chaplygin gas [9] 
and tachyon models [10], dark energy and dark matter 
are described as separate limiting cases from an exist- 
ing unique field. Partly inspired by these models one 
could naively assume the existence of a unique axion field 
which, when excited, would make dark matter, and when 
at its vacuum ground state would be the source of phan- 
tom energy. Coincidence time could then be interpreted 
as the time when both the vacuum and the excited states 
are approximately equally populated. Of course, this 
would not solve the coincidence and fine tuning problems 
but provided some explanation to these problems and to 
the prediction of generating such strangely small amounts 
of homogeneously distributed axions. It remains never- 
theless an intriguing possibility that such a small value 
of the cosmological constant be supplied by the potential 
energy density of an " ultrainvisible" axion field which 
can be dubbed quintessence axion [22]. 
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III. THE COSMOLOGICAL MODEL 



IV. THE FIELD THEORY 



We shall consider in what follows a simple model 
where the massless scalar field is a quintessential field, 
equipped with a potential V(<fi), which is minimally cou- 
pled to Hilbert-Einstcin gravity. The action for this 
model is 
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where L m again is the Lagrangian for observable matter 
fields. If we again restrict ourselves to the case where 
(i) L m — 0, (ii) both the scalar field and Hilbert-Einstein 
gravity satisfy the Fricdmann-Robertson- Walker symme- 
try, and (iii) we assume a perfect fluid equation of state, 
then by integrating the conservation law for cosmic en- 
ergy, dp = — 3(p + p)da/a, it turns out that the energy 
density for the quintessence scalar field will be given by 



(8) 



in which R is an integration constant. From the Fried- 
mann equation derived from action (7) subject to the 
FRW symmetry and L m — 0, 
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where ao and to are the initial radius and time, respec- 
tively. Note that for ui > —1 this solution describes 
an accelerating universe whose scale factor increases to- 
wards infinity as t — > oo. We are here most interested in 
the observationally not excluded yet case where u> < — 1 
which corresponds to the so-called phantom dark energy 
for which the dominant energy condition is generally vi- 
olated, i.e. [15] 



p + p<0, 



(11) 



even though the energy density is surprisingly ever in- 
creasing. Notice furthermore that in this case the scale 
factor blows up at a finite time, 
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giving rise to a "big rip" [15]. One would expect that 
in the vecinity of time £* the size of all local geometrical 
objects increased in a rather dramatic way. We see that 
the larger ao and \ui\ the nearer the doomsday, and hence 
the time expected to have macroscopically increased ge- 
ometrical objects in our universe. 



Now we shall consider the field theory which is asso- 
ciated with phantom energy, starting with the general 
formalism for any value of lo. From the equation of mo- 
tion for the scalar field 6, 
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where H = a/a and V((f>) again is the potential for the 
field (f>, and the expression for kinetic energy of the field 
[Eq. (2)], we can obtain a relation between dV/d(p and 
the scale factor which can be expressed in the form 



■ t y/A(l+u)(l - u)a 



-3(l+w) 



dV 



(14) 



On the other hand, integrating the equation of motion 
for the field (f> we can derive the dependence of this field 
with time t. 
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where l v = y/R/A is the Planck length and B is an 
integration constant. 

From Eqs. (10) and (15) an expression for the scale 
factor in terms of the scalar field can now be derived. It 



is: 
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Using Eqs. (14) and (16) we finally obtain the expression 
of the field potential for any ui < —1/3 as a function of 
the field 6 



v{4) = v +±(i- 



-w)VfiB ~ 2 exp (-3V1 + uj£ P 4>) , (17) 



in which V is a constant. Potentials with this form nat- 
urally arise in supergravity models [23] and have been 
used in a variety of context, ranging from accelerating 
expansion models [24] to cosmological scaling solutions 
[25]. 

Together with the expression for the field <j> in terms 
of time t, this expression solves the field-theory problem 
for any value of u) < —1/3. According to our discussion 
above, for axionic phantom energy we have to perform 
the continuation (f> — > i$ in addition to taking w < — 1. 
In such a case, the solution of the field theory problem 
would read 

*(*) = 
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FIG. 1: The potential V for imaginary (solid lines) and real 
(dashed line) scalar field. The scales for the potential and time 
are arbitrary, t* is the time at which the big-rip singularity 
would occur in the future. It is worth noticing that the singu- 
larity in the spacetime curvature at t* coincides with a singu- 
larity in the imaginary field <E> and the potential V($). In con- 
structing this figure it has been assumed that the initial scale 
factor and the parameter entering the respective equations of 
state are such that a » [(I - w d )/(l + |w p |)] 1/(3(wd+|wpl)1 , 
where ui d > — 1 and lo p < —I, with ao sufficiently large. 
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(19) 

As to the physical motivation for the exponential ax- 
ion potential (19), I will briefly comment on the influence 
that the matter-field sources may have on the stability 
of the present axion model of phantom energy If, corre- 
sponding to the physically interesting case that V m ^ 0, 
we would let the influence of matter sources on the expan- 
sion factor immediately before of phantom domination to 
have the form of a perturbation, and analyze the result- 
ing model via phase space [26] , it can be seen that if we 
choose Vq = then potential (19) and the solution given 
by Eq. (10) (with w < -1) and Eq. (18), that is for 
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for phantom energy domination and catastrophic big rip 
with u> < — 1 and 
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for which oj = — 1 — A^/3 (See also phantom models with 
Born- Infold type Lagrangians [27].) 

Note that the potential V(<I>) takes on a value 



V = V + hi + VRa^-V 



(20) 



at the initial time t = t to steadily increase thereafter up 
to infinity when t = t* , at the big rip. From that moment 
on the potential would continuously decrease, tending to 
reach its constant minimum value Vo as t — > oo (see solid 
curves in Fig. 1). That behaviour is in sharp contrast 
with that of potential V(<j>) for dark energy with w > — 1, 
which starts with a value 



V = V 



)VRa a 



3(l+u;) 



to monotonously decrease down to the value Vo as t — > oo, 
such as it is also shown by the dashed curve in Fig. 1. 
We have solved in this way the field theory which is as- 
sociated with dark and phantom energy if the latter is 
interpreted as originating from the existence of an ax- 
ionic field which is expressible as a massless, purely imag- 
inary scalar field. The explicit emergence of axions in the 
cosmological context of phantom dark energy appears to 
be just another example where such rather elusive par- 
ticles are invoked as candidates for inexorably needed 
cosmic constituents which are themselves defined to be 
unobscrvablc. Actually, apart of being a key ingredi- 
ent to solve the strong-CP problem [20] or to represent 
quantum-gravity topology changes [19], axions have al- 
ready been widely used to solve a variety of astrophysical 
and cosmological shortcomings [28]. 

It appears of some interest to briefly comment next 
on a potentially fruitful consequence from the connec- 
tion of axions with our cosmological model. If we are 
actually living in a universe with phantom energy where 
there will occur a big rip singularity which can never be 
circumvented by causally- violating connections to the fu- 
ture, then the contracting branch of solution (9) can be 
physically disregarded as the big-rip singularity is a true 
curvature singularity. In such a case, only the branch 
of potential V($) before the singularity can have phys- 
ical significance, and hence the global minimum of that 
potential would become placed at t = to with the value 
given by Eq. (20). It corresponds to a value of the field 



correspond (see e.g. Ref. [26]) to a dynamical attractor 
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Interpreting potential (19) as that axion potential result- 
ing from the QCD nonperturbative effects, we can take 
[29] <f> m = Ja^cS, with Ja a constant and 9 e s the ef- 
fective 6 parameter resulting after the diagonalization 
of the quark masses in the QCD Lagrangian. Setting 

then f A = 2/[3£ p ^J\^l], 9 cS = In ^B al iH - 1)/2 ^ and 

Bq = 1 it follows that, since 6 c s < 10~ 9 [29], the initial 
radius of the phantom energy universe ao should be very 
close to unity, that is the onset of such a phantom energy 
regime, if it ever at all occurred, must always be placed 
at nearly just the present epoch. Actually, if 9 c ff is can- 
celled to completely solve the strong CP problem [29], 
then the choice Bq = 1 implied that ao would exactly 
satisfy ao = 1. In such a case, the observational setting 
of the value of the state equation parameter w would also 
set the value of ao, and hence of time t*. 

A potential problem with the proposal in this paper 
could be that whereas current axion theories have a cut 
off on the scale of inflationary energy, typically at around 
the GUT characteristic energy of 10 16 GeV [29], in order 
to avoid phantom energy decaying into gravitons, a cut 
off of < 100 MeV should be introduced in the effective 
phantom theory [30]. Not with standing, while the for- 
mer cut off appears relevant for dark matter axions which 
can initially be highly excited, the axion field proposed 
in this paper as the source of phantom energy can never 
be excited outside the vacuum. On the other hand, if we 
keep the current values for color anomaly of the Peccei- 
Quinn symmetry, the pion and quark masses and the pion 
decay constant, then the above Peccei-Quinn symmetry- 



breaking scale, Ja ~ 10 19 (|w| — 1) 1 GeV, would imply 
a mass of the phantom axions 



m A ~ 0.62 eV " ~ lO^/M" 1 eV, 

JA 

i.e. just at the extremal minimum of the allowed values 
for the axion mass, and clearly well below the cut off 
required by Carroll et al. [30] for phantom energy to be 
stable. 



V. SUMMARY 

Imposing the weak energy condition, the negative ki- 
netic energy for a phantom field is interpreted as being 
originated from super-light axions. We have built up a 
simple cosmological model encompassing constant equa- 
tions of state for both oj > — 1 and u> < — 1. In the latter 
case the universe evolves toward a singularity at finite 
time. The field theory is then solved for all cases, deriv- 
ing an increasing exponential potential for the phantom 
field which also diverges at the singularity and shows a 
dynamical attractor also when matter fields are present 
that inexorably leads to a catastrophic big rip. 
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